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Abstract 

For thermal systems, standard perturbation theory breaks down be- 
cause of the absence of stable, observable asymptotic states. We show, how 
the introduction of statistical quasi-particles (stable, but not observable) 
gives rise to a consistent description. Statistical and spectral information 
can be cleanly separated also for interacting systems. 
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1 Introduction 



It is a well established fact by now, that naive perturbation theory breaks down 
at finite temperature, see refs. ||l|, |^ for a recent overview. This is due the 
modification of space-time symmetry in the presence of matter or a heat bath 

i.e. to the absence of stable asymptotic states for the observable physical 
particles. To overcome this problem, first of all one has to employ a "doubling" 
of the Hilbert space, leading to two-point Green functions that are 2x2 matrices. 

The Schwinger-Keldysh jj] or Closed-Time-Path method (CTP) is such a 
model, which has been used to derive transport equations, e.g. for nuclear 
phenomena, over several decades ||^, ^. However, the physical interpretation of 
the matrix structure of the propagator remains quite obscure in this formalism, 
and no justification is obtained for a perturbation expansion. 

Another model is called Thermo Field Dynamics (TFD) §, |. It differs 
from CTP in the important aspect, that the doubling of the Hilbert space is 
due to two disjoint, mutually (anti-) commuting representations of the canonical 
commutation relations. Thus, in contrast to CTP, TFD has a firm mathematical 
basis and can be used perturbatively. To this end, one has to sacrifice either 
the stability aspect or the observability aspect of asymptotic states. The first 
leads to a perturbation expansion in terms of unstable particles [^] . It is the 
purpose of the present paper to demonstrate, that the second way leads to a 
simpler and straightforward approach. 

To this end we show, that the thermal instability of observable states can be 
absorbed into a Bogoliubov transformation also for interacting systems. This 
Bogoliubov transformation then defines stable, albeit non-observable, quasi- 
particles which serve as basis for a perturbation expansion Throughout 
this paper we term them, for historical more than physical reasons, statistical 
quasi-particles |T^. 

2 Thermal Bogoliubov transformation 

To establish the notation, we first discuss TFD for a single bosonic quantum 
state. We introduce creation and annihilation operators a\ a, a\ a for the two 
representations, with canonical commutation relations 

[a,at]=l [a,a^=l [a,a^^O (1) 
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(see ref. [g| for a complete discussion). The a, a operators annihilate the 
physical vacuum ||0, 0))and the two sets are transformed into another by means 
of an anti-unitary mapping, called the tilde conjugation (see [p^ for the tilde 
conjugation rules). For brevity, we work in the a = 1 representation, i.e., the 
thermal equilibrium state at inverse temperature /3 and chemical potential /i is 
described by two state vectors 

worn = exp(/atSt) ||o,0)) / ^ g-'^^-'^) 

(2) 

{{Om\ = ((0,0|| exp(aS) 



Within this framework, the ensemble average of an observable £[a, a^ is calcu- 
lated as the expectation value 

\ ^ ({om\£\\om 

I iomom ■ ^ ^ 

Obviously then the state vectors ((0(/3)|| and ||0(/3))) are annihilated by certain 
linear combinations of the above operators, 

mm^^ eiio(/3))) = o ((o(/3)iiett = o io{m'^ = ^. (4) 

obtained as 



(5) 



IB is a 2x2 matrix with determinant 1. Since the ^-operators obey similar 
commutation relations as do the a-operators, they define our quasi-particles. 
It is crucial to realize, that these entities have nothing in common with the 
"usual" definition of quasi-particles, which refers to physical states with an 
almost pointlike mass spectrum. 

Eqn. (^) is essentially a Bogoliubov transformation |l4j. The most general 
form for the Bogoliubov matrix compatible with our choice of state vectors is 
then 

B.-^„p(,.3)(_\ -/) , (6) 
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where T3 = diag(l, —1) is a Pauli matrix, s is free parameter and / is the sta- 
tistical weight of a physical single particle state within the ensemble. Therefore 
the number density of the physical particles is 

" (7) 



1-/ 



We find, that the following considerations are extremely simplified by chosing 
s — 1/2 log(l + n) = —1/2 log(l — /), since then the Bogoliubov transformation 
is linear in the density parameter 




3 Interacting Boson field 

We now consider a fully interacting bosonic quantum field. At finite temperature 
the irreducible representations of the space-time symmetry group are character- 
ized by two rather than one continuous parameter Hence, the interacting 
field can be expand into modes with definite energy and momentum as 

oc 

^{x)= JdE J^=p'^^iE,k)<PEMi^) (9) 

The commutation relation of these fields are, in general, not known. However, 
we want to calculate only the two-point Green function of the interacting field, 
i.e. it is sufficient to know the expectation value of the commutator of these 
fields. 

This expectation value in turn can be absorbed in the definition of the weight 
function p{E,'k), i.e. for the field operators we can define 

,E.M{t.^),dt^E',i.'{t,^')]) = e'^'^^-^'HESiE - E')6{k-k') . (10) 



In other words, for the calculation of bilinear expectation values of interacting 
fields it is sufficient to consider the (^_B,k(a;) as generalized free fields Q. The 
full information about the single-particle spectrum of the theory is contained in 
the weight function p{E,'k), and we require the normalization 

JdE^ p{E,k) = 1 . (11) 
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Note, that the existence of this spectral decomposition is only guaranteed in 
case the system is space-time translation invariant, i.e., if it is in a thermal 
equilibrium state. 

We now apply the thermal quasi-particle concept to each energy-momentum 
eigenmode of the system, i.e. we define quasi-particle operators S,E,k associated 
with energy-momentum eigenstates and annihilating the statistical state vectors 
as in (^). These are then Bogoliubov transformed into physical particle operators 
for definite energy and momentum, and those arc summed with the above weight 
function to operators ak{t) such that the interacting field is 

cl^ix) = J {al{t)e~''^- + a,it)e''^-) (12) 

(note the absorption of the usual energy normalization factor into the operators). 
In this expansion of course, the "creation" and "annihilation" operators are quite 
complex objects. From the above reasoning, i.e. their decomposition into the 
modes of definite energy and momentum, we then obtain 

{{Om\al^''\t) 4^)(t')||0(/3))) 

= 53(k_k') JdEpiE,k) {T,B^inEi.)Ts)^^{B^inEu));,'e^^^'-''^ 

= 6^k-k') jdEp{E,k) (S(nsk)):i'(r3S(nsk)T3)i,e-i^(*-*') . 

(13) 

The parameter n^jk defining the individual Bogoliubov transformations thus ap- 
pears under the energy integral. While this result can be understood intuitively, 
i.e. every mode is in thermal equilibrium with the rest of the system, it can also 
be understood in a more formal way. 

To this end on has to look at the time evolution of initially free particle 
operators: it is highly nonlinear. Therefore the Bogoliubov transformation in 
interacting systems is non-linear, and this non-linearity is reflected in the above 
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energy integral. It is indeed possible to derive the above equation without ever 
touching the concept of generalized free fields . 

Note, that in our notation the weight function p has only support for positive 
energy arguments. The retarded and advanced propagator are in momentum 
space 

i:)^"4(£;,k) = fdE' p(E',k) ( ^ , (14) 

^ ' ' J ' ' \E - E' ±ie E + E'±ie) ' ^ ' 

and the limit of free particles is recovered when 

piE,k)^SiE'-u;l)eiE) . (15) 
We then obtain for the propagator matrix 



D^''''\t,t';k) = -i JdEp{E,k) 



+ r^B^in^^) (^^^''^'^ ^e(i-o) (^^("^''»"^'''^"*'') • 

(16) 

This is a straightforward generalization of the result from ref . ||l^ . In the free- 
particle limit, the textbook result for the finite temperature boson propagator 
is recovered. 

The above expression is still somewhat unsatisfactory, since the i3-matrices 
are subject to an energy integration. However, because of the special form (||) 
we chose for the parametrization, the integrand is linear in the parameter n^k- 
Thus the integration can be carried out if one defines 



= jdEp[E,k)nE,. (e-i^(*-*')+e 

Z{t,t') = ydi?p(i?,k) (e-i^(*-*')+ei^(*-*')) 



\E(t-t'] 



(17) 
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Some elementary matrix operations then lead to the result for the propapagator 
00 



= -iZ{t,t'){B{N{t,t')))-^ (^ii-f) B{N{t,t'))T, 

(18) 

Here we have kept positive and negative energy states separate: the inner prop- 
agator matrices are diagonal in this case. One can, however, also combine the 
two parts into one triangular inner matrix, sandwiched between two Bogoliubov 
matrices: 



{B{N{t,t')))-^ X 

(19) 

(-iQ{t - t') {Z{t, t') - Z*{t, t')) i (1 + 2N{t, t')) Z*{t, t')\ 

\ ieit'-t)izit,t')-z*it,t'))) 

X B{N{t,t'))T3 . 



The physical relevance of the function N{t,t') diagonalizing the propagator 
becomes obvious, when we consider its equal time limit. It approaches a constant 
then, 

l,hn N{t, t') = iVf Ihn ^7V(t, <') - . (20) 
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Comparision to ( p^ gives 



ldEp{E,k)nEi. lim / d^k ((0(/3)||4 (t)afe, (i')||0(/3))) 

ivf = = 7 • (21) 

dEp{E,k) dEp{E,k) 

In other terms, is the time- independent equihbrium Heisenberg density of 
the physical particles with momentum k. The quantity N(t,t') therefore is the 
observable fluctuating particle number of these modes. 

The separate diagonality of the inner matrices of the propagator, i.e. the 
requirement of unperturbed statistical quasi-particle propagation, is therefore 
equivalent to chosing the correct physical Bogoliubov parameter. 



4 Conclusions 

The concept of quasi-particles in statistical physics is known for some time , 
and it is also known that a linear relation between the matrix elements of an 
interacting propagator can be used to bring it to a triangular 2x2 matrix form 
[ p!?] . We have put these two things together and introduced statistical quasi- 
particles into Thermo Field Dynamics. This leads to a diagonal propagator for 
nonrelativistic models Q. 

When negative energy states are taken into account, the full propagator can 
also be written as diagonal matrix sandwiched among Bogoliubov matrices, but 
separately so for particle and anti-particle states Their combination then 
gives a triangular inner propagator matrix, see eqn. (|l9|). We find, that the 
Bogoliubov transformation necessary for this diagonalization is given in terms 
of a physical parameter, the observable fluctuating particle density N{t, t'). This 
diagonalization also defines stable physical modes with fixed momentum and an 
average energy Q . 

While this is clearly a conceptual advantage, the diagonalization of the full 
propagator at finite temperature also has a tremendous technical advantage over 
the CTP formalism. In effect our method separates the statistical information 
about the system (boundary conditions, thermal particle-hole excitation etc.) 
from the purely spectral information contained in the weight function k). 
Let us note, that the same separation can be achieved for non-equilibrium sys- 
tems. 
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The application of the quasi-particle concept to a given system not only 
demonstrates this technical advantage, but lends a physical interpretation of 
the 2x2 matrix structure of thermal quantum theories |Q, |^. For brevity we 
only state the results: Requiring, that the triangular propagator (|l8|) solves the 
diagonal components of a Schwinger-Dyson equation gives p{E, k) as function 
of real and imaginary part of a retarded self energy function. 

The off-diagonal component of the Schwinger-Dyson equation contains the 
statistical information, i.e. it is a consistency criterion for the function N{t, t'). 
For the timc-indcpcndcnt case considered here, we were able to derive this 
consistency criterion as the condition of global thermal equilibrium. For non- 
equilibrium systems, where a similar separation of statistical and "spectral" 
information can be obtained, the diagonalization condition for the propagator 



is nothing but a transport equation 1 1 1 
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